Quantum phase transitions in photonic cavities with two-level systems 
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Systems of coupled photonic cavities have been predicted to exhibit quantum phase transitions 
by analogy with the Hubbard model. To this end, we have studied topologies of a few (up to six) 
photonic cavities each containing a single two-level system. Quantum phase space diagrams are 
produced for these systems, and compared to mean-field results. We also consider finite effective 
temperature, and compare this to the notion of disorder. We find the extent of the Mott lobes 
shrink analogously to the conventional Bose-Hubbard model. 



I. INTRODUCTION 

There has recently been a convergence of several dif- 
ferent fields of physics: condensed matter, quantum op- 
tics and information science. This convergence has been 
realized by a staggering increase in the ability to fabri- 
cate and control quantum systems experimentally, and 
an ability to attack theoretical problems of increasing 
complexity. One aspect of this convergence of fields is 
the push to realize the quantum computer. Here we 
discuss another aspect: that of a quantum simulator. 
In particular, we explore the possibility for a quantum 
atomic-optical system (here an interacting lattice of op- 
tical cavities with embedded two-state systems) to un- 
dergo a quantum phase transition by direct analogy with 
the Hubbard model. 

The Hubbard model [l[ describes the hopping of in- 
teracting particles around a lattice of allowed positional 
states. A quantum phase transition is observed be- 
tween delocalized particles (supcrfluid phase) and local- 
ized particles (Mott-insulator phase) depending on the 
strength of the hopping term relative to the onsite in- 
teraction. Numerous facets of the Hubbard model have 
been considered including the prediction of glassy phases 
0, Hilbert-space optimization d, and implementations 
of topological quantum computing 0] . One of the most 
dramatic and beautiful examples of the Bosc-Hubbard 
model is the prediction [f| and recent demonstration Q 
of the quantum phase transition in an ultra-cold atomic 
gas. Such demonstrations are significant for applying 
canonical solid-state treatments to the more controllable 
regime of atom optics, allowing new predictions to be 
tested (e.g. the supersolid phase 0]). Recent work on 
quantum phase transitions in photonic band-gap lattices 
does the same for the photonics-solid-state boundary 

I IE, El El El. 

At first glance, the possibility for a quantum phase 
transition in an optical system seems surprising. This 
is due to the fact that photons do not normally inter- 
act with each other with any appreciable strength, and 
for this reason most nonlinear optical processes are con- 
fined to the realm of classical optics. There are, how- 
ever, many exceptions to this, but perhaps the most dra- 



matic is the phenomenon of photon blockade [l4|. Pho- 
ton blockade is an example of a cavity quantum electro- 
dynamical interaction in the strong coupling limit. An 
atom is placed in a cavity, and because the energy lev- 
els of the atom-cavity system depend on the number of 
photons in the cavity, aphoton-number dependent reso- 
nance shift is observed [l5[ . If the atom-cavity interaction 
is strong enough, this shift can be sufficient to prevent 
more than a predetermined number of photons to enter 
the cavity: photon blockade. This effect has been an- 
alyzed for four-state systems [3, EH, E3 an( l two-state 
systems 1181. More recently photon-blockade has been 
observed [l9l [2fj| . adding substantial impetus to apply 
this effect to a range of applications. 

Here we consider the properties of a lattice of cavities, 
each containing a single, quasi-resonant two-state system, 
so as to be effectively treated by the Jaynes-Cummings 
interaction [2l[ . We go beyond the earlier idealized treat- 
ment @ by building systems of increasing size to predict 
the results in the few-cavity (up to six) limit, and also 
consider the thermodynamic implications of disorder. By 
directly connecting these small scale cases (solved by di- 
rect diagonalization) with the thermodynamic limit, our 
results serve as a guide to coming proof-of-concept ex- 
periments. We are also able to compare our finite cases 
with the thermodynamic limit @, El S3, H, EEL HI ■ 

In Sec. [IT] we introduce the system of coupled photonic 
blockade cavities that will be investigated for quantum 
phase transitions, and the extended Jaynes-Cummings- 
Hubbard Hamiltonians for both the exact calculation and 
for the mean-field approximation. In Sec. IIIII we present 
results from the exact diagonalization techniques and 
compare with mean-field solutions. Finally, we consider 
disorder and the implications for effective model temper- 
ature in Sec. IIVI 



II. PHASE TRANSITIONS IN THE 
JAYNES-CUMMINGS-HUBBARD MODEL 

The system under consideration is a lattice of optical 
cavities, each containing a single quasi-resonant two-state 
system. The canonical treatment for a single atom-cavity 
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FIG. 1: (Color online) This schematic shows a possible im- 
plementation of the system, for three nearest neighbors. The 
dielectric medium is shown in yellow. The photonic crystal 
is made by periodic variations in refractive index, caused by 
drilling holes (white disks). The cavities are regions where 
holes have not been drilled - effectively, the undrilled holes. 
In each cavity, a red disk represents the two- level atom. The 
arrows indicate the nearest neighbors, across which photons 
can tunnel with hopping rate k. 



system is the Jaynes-Cummings model. Photon hopping 
between cavities (which is effected by leakage out of the 
cavities, and into neighboring cavities) allows the direct 
comparison to Hubbard systems p7j . and hence we re- 
fer to this as a Jaynes-Cummings-Hubbard (JCH) model. 
There are numerous ways in which to realize such a sys- 
tem depending on the available experimental configura- 
tions and desired topologies; for example, photonic band 
gap structures @ and coupled-cavity waveguides 
perhaps realized in micro-fabricated diamond [2a 



arrays of superconducting strip-line cavities [3ll . [32 
microcavities with individual cold-atoms connected via 
optical fiber interconnects [HI, or plasmonics [HJ|. For 
concreteness, we will focus our attention on a photonic 
bandgap structure where a two-dimensional array of pho- 
tonic bandgap cavities constitutes the underlying lattice 
and defines the nearest-neighbor topology, and the two- 
state system is realized by an implanted impurity, Fig. [1] 

Our emphasis in this work is on systems with a sin- 
gle two-state system per cavity, but it is important to 
note that this is not the only potential system for ob- 
serving similar quantum phase transitions. Hartmann et 
al. have considered four-state systems [a. [35l| . in keeping 
with the original Imamoglu proposal [l4j . whereas the 
case of many atoms per cavity has been considered by 
Na et al. [ll[ . The approach of Ref . [ll[ is particularly 
useful for providing a clear path to experiments using 
GaAs quantum dots. 

To understand the properties of the JCH system, we 
first review the properties of the individual atom-cavity 



(Jaynes-Cummings) interactions. The Hamiltonian is 

W JC = e(7+cr~ + cja f a + j3(a + a + o-~a f ), (1) 

where er + and o~~ (a^ and a) correspond to the atomic 
(photonic) raising and lowering operators, respectively. 
The transition energy of the atomic system is e, the cav- 
ity resonance is u>, and the cavity-mediated atom-photon 
coupling is (3. The difference A = to — e is the detuning. 

Let \g, n) (|e, n)) (n G Z*) represent a cavity that con- 
tains n photons and a single two-level atom in the ground 
(excited) state. The energy eigenvectors of Eq. (Q]) arc 
given by \g, 0) and 



\±,n) 



with eigenvalues 



-(A/2)± X (n) 



1) 



v/2x 2 (n)Tx(«)A 



(2) 



Vn > 1, 



E\ g ,o) = 0, £|±,„) = nuj ± X (n) - A/2, (3) 
where we have used the generalized Rabi frequency 



X(n) = + A 2 /4 Vn > 1 



(4) 



where n is the total number of excitations. These eigen- 
states correspond to the well known dressed (polaritonic) 
states, and we call the basis formed by them the single- 
cavity dressed basis. The eigenspectrum of a single atom- 
cavity system is shown in Fig. O Because of the atom- 
photon induced shift of the energy levels as a function 
of the number of excitations in the system, there is an 
effective photon-photon repulsion [l^]. It is this photon- 
photon repulsion which plays the role of the on-site term 
in the Hubbard model, however it is important to note 
that because the repulsion decreases with an increasing 
number of particles, the canonical Bosc-Hubbard system 
is not realized in our case, and so although many qual- 
itative similarities are predicted between the JCH and 
Bose-Hubbard models, exact equivalence is not guaran- 
teed. 

The non-bosonic nature of the particles in the Hamil- 
tonian of Eq. |T]) requires further discussion. Neither the 
JCH system nor the four state system with few atoms 
per cavity 0] retrieve bosonic commutation relations. 
The limits where we can view the system as being com- 
prised of interacting bosons arc many atoms per cavity 
(holds for both the Jaynes-Cummings and four state sys- 
tems BtniUll), large det uning [9(, and large excitation 
number. Arguably the most important case is that de- 
scribed here, namely one atom per cavity with few exci- 
tations. This is because this regime maximizes the non- 
linear (photon-photon) interactions, and is therefore the 
most experimentally accessible regime. 

Differences between the JCH and Bose Hubbard sys- 
tems are interesting topics for investigation, and a study 
of the particle nature should prove fruitful, but goes 
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L = ~}2d = i Li, where Li = cr^cr^ + a\di is the number 
operator of atomic and photonic excitations of the zth 
cavity. One can include a term —fiL in the Hamilto- 
nian. We show below, through arguments of statistical 
mechanics, that \x represents the chemical potential. Sec- 
tion [IV] will continue into a discussion on effective model 
temperature and disorder. Let us include this chemical 
potential term in the Hamiltonian directly as follows: 
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FIG. 2: (Color online) Energy eigenvalues [Eq. Q] for a single 
cavity as a function of the detuning A. We set uj — 12/3 
for illustrative purposes. The blue dashed lines indicate the 
asymptotes of the energies for each band, the diagonal lines 
represent the energy of the cavity due to the atom, while the 
horizontal lines represent the energy due to the photons only. 



beyond our present work. We may understand some 
of the differences by comparing the onsite repulsion in 
the Bose Hubbard and JCH cases. In the Bose Hub- 
bard system, the interaction U is a constant, however 
in the JCH model this can be seen as having a parti- 
cle number dependence, i.e. U±(n) = x( n + 1) — x( n )- 
In the large photon limit, we obtain a non- interacting 
Bose g U±(n) — > 0, and in the large detuning limit, 
U±(n) — > ±/3 2 /A, which is a constant bosonic Hubbard 
type repulsion [TJ HB] ■ There is also no ideal Kerr- type 
term to generate an exact quartic interaction. Nonethe- 
less, as has been shown, qualitative similarity between 
the phase diagrams of JCH and Bosc-Hubbard systems 
is found, and the analysis of these phase diagrams is a 
major topic of this paper. 

To generate the JCH Hamiltonian, we add hopping 
between cavities, and for a system of N cavities we have 



N 



(id) 



(5) 



where individual Jaynes-Cummings Hamiltonians of 
Eq. ([1]) have identical e, u> and (3, (this restriction will 
be relaxed later). The intercavity hopping occurs with 
frequency Kij = k for nearest neighbors, and k = other- 
wise, it is this term which defines the topology of the net- 
work. Photon transmission through a one-dimensional 
chain in a similar structure has also been considered [3?) • 
To divine the properties of the phase transition seen in 
the thermodynamic limit, we introduce the operator that 
measures the total number of excitations of the system 



n' = n- nL. 



(6) 



We assume that the entire N cavity system with I total 
excitations exists in the ground energy cigenstate |?/> g ), 
so that W'lV'g) = Egli^g) an d L\ip g ) = l\if) g ) (i.e. these 
two operators commute). 

To show that fi has the general form of a chemical po- 
tential, we begin with the usual definition of free energy 
F = E — TS, where E is the energy of the system Ti (be- 
fore chemical potential has been included), T is temper- 
ature and S is entropy. Assuming that T = and using 
the definition of chemical potential as the derivative of 
the free energy with respect to number of excitations 



dF 
~dl 



T,V 



dE 
~dl 



(7) 



T,V 



where / is used, as excitations act like particles in this 
system. We use the Hellmann-Feynman theorem [HI, |39| 
to calculate the derivative of the energy with respect to 
number of excitations 



dE an 

= (Vg||(W' + M0lV'g), 
= M- 



(8) 



Hence the [i of Eq. ^ represents chemical potential as 
required. 

The TV cavity bare basis consists of state vectors of 
the form |si,m) ® |s 2 , n 2 ) ® • ■ ■ ® \s N , n N ), s l G {g, e}, 
rii G {0, 1, . ..}. In principle, this basis is infinite in ex- 
tent, because the number of photonic excitations per cav- 
ity is unbounded. By ordering the bare basis by the total 
number of excitations (either photonic or atomic) across 
all cavities, one may express the Hamiltonian of Eq. ([5]) 
in block diagonal form H = diag[H (0) ,H (1) ,H (2) , . . .], 
where T~ft) is the matrix corresponding to / excitations. 
The size of each block is determined by the number of 
ways in which the excitations can be shared between the 
atomic and photonic degrees of freedom. We denote the 
number of states (equal to the size of the matrix TL^) as 
s, where 



min(Z,JV) 

E 



s. 



N 



(9) 
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The above summation has two terms. ( .) is the total 
number of atomic excitations across the lattice (note that 
on each site the number of atomic excitations can only be 
zero or one) . represents the number of photonic exci- 
tations, and is the number of ways to share the I — i pho- 
tons between the N cavities, (e.g. 5| = length[(2, 0, 0), 
(0,2,0), (0,0,2), (1,1,0), (1,0,1), (0,1,1)] = 6). 

To gain insight into the problem, we explicitly show 
and TiW for a two-cavity system in the bare basis 
as follows: 



H 



(o) 







H 



(i) _ 



oj (3 



-K 





(10) 



-k w 13 
(3 e 



where the (two-cavity) basis for Tt^ is {\g,0, g,0)} and 
for is {\g,l,g,0),\e,0,g,0),\g,0,g,l),\g,0,e,0)}. _ 

To connect the finite case with the thermodynamic 
limit, we examine the phases of the N cavity system 7i' 
[Eq. ijBJ)] as a function of u>, k, and A. In particular, 
we are concerned with the expectation value of the total 
number of excitations of the system (L) . Note that L is 
diagonal when represented in either the dressed or bare 
basis (but has different values in each) . From , by sub- 
tracting fiL from our Hamiltonian Eq. (EJ), we determine 
that (L) = —dE g /dfi, where E g is the ground state en- 
ergy of the extended Hamiltonian TC . Some preliminary 
analytics can simplify the calculation of (L) considerably; 
we show this now. 

We begin by noting that the part of L corresponding 
to exactly I excitations, represented as ZA", has the very 
simple form 



L® = ll. (11) 
Again employing the Hcllmann-Feynman theorem, 



dE, 



d 



-<^l^(« (0 -AdlM>, (12) 

where E^ is the ground state energy, and is the 
corresponding eigenstate, of TL^ — fiL"'. This reduces 
to 



dE, 



(/) 



= -I. 



(13) 



So, if 



M = { min[cigenvalues(H 



(0) 



L(°))], 



min [eigenvalues (7i { - L ' 1 — /iL*- 1 -*)], . . .} 



(14) 
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FIG. 3: (Color online) Ground state energies for each block 
of the Hamiltonian, Eq. as a function of (fx — uj)//3, with 
k = 0. Recall that (L) is equal to the negative of the slope 
(with respect to /i) of the smallest energy eigenvalue of the 
whole Hamiltonian, and that this slope is always an integer. 
This figure therefore shows explicitly how the phase diagrams 
are constructed - the black points indicate the boundaries 
between plateaus. 



and 



then 



/: {0. 1....!' -M. 



(L) = /- 1 [min(M)]. 



(15) 



(16) 



In short, to find (L), one simply needs to locate which 
block the minimum eigenvalue of [TL — fiL) corresponds 
to. Obviously, (L) can only have non-negative integer 
values. This is illustrated further in Fig. [3] In this figure, 
the smallest eigenvalue of ,H^\ ... is plotted as a 
function of (p, — u))/P for n/[3 = 0. For each value, (L) is 
given by the negative slope of the smallest eigenvalue at 
that point. 

We now introduce the mean-field Hamiltonian. The 
mean-field approximation focuses attention on one par- 
ticular cavity, and assumes that its z nearest neighbors 
(that is, the coordination number is z) all behave like it. 
To invoke the mean-field, we use the decoupling approx- 
imation a\aj = (a\)aj + (aj)a\ — (a\)(aj), and introduce 
the superfluid order parameter ip = (cti) = (aj) (which 
we assume real), so that the Hamiltonian of Eqs. ([5]) and 
(El) becomes 



H 



MF 



H 



JC 



ZKip((r + a) + zn'ij} 2 — /iL. 



(17) 
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The basis uses just one cavity, but the system (approxi- 
mately) describes an infinite number. Note that the num- 
ber of nearest neighbors z effectively "renormalizes" the 
mean- field coupling, i.e., K — > zk. 

We have also considered using a larger unit cell - for ex- 
ample, using two cavities, each with z nearest neighbors. 
We found that while this is more difficult to calculate 
(finding eigenvalues of larger matrices) this exactly repli- 
cates the results of the original mean-field - which is not 
surprising. However, this technique could be used to in- 
clude disorder in the infinite cavity limit. Note that while 
Eq. (|16p informs that the total number of excitations is 
integer, the mean-field Hamiltonian of Eq. (|17[) informs 
of the number of excitations per cavity. Accordingly, we 
will find equivalence when the number of excitations is a 
multiple of the number of cavities. 



III. RESULTS 

In this section we analyze the quantum phase diagram 
of the Hamiltonian of Eq. §5§ for various topologies. Our 
analysis is based on exact diagonalization for up to six 
cavities. Topology is implemented through the Kij terms 
of Eq. ([5]). We then compare these topologies with the 
mean-field approximation. 

We begin by considering the quantum phase diagrams 
of the exact systems. We display the phase diagrams 
corresponding to two, three, four, and five cavities ar- 
ranged in one dimension with periodic boundary condi- 
tions in Fig. 0J Each color corresponds to a different 
plateau, a constant state in excitation space - these are 
Mott insulating phases. It is worth pointing out that in 
our discrete model, no supcrfluid phase exists. However, 
for significantly large k, the plateaus get closer together, 
approximating the superfluid phase diagram, as in the 
mean-field case. 

In total, eleven topologies were examined. These are 
listed in the first column of Table |TJ The topologies of 
a square, triangle and six cavities with z = 3 could be 
considered special, as they can represent infinite square, 
triangular and hexagonal lattices respectively. However, 
no significant differences (with respect to matching of 
phase diagrams to mean-field) are found between these 
topologies and the rest. 

For all topologies, a "pinch" effect is noted as k — > 
between (L) = N and (L) = 2N, between (L) = 2N 
and (L) = 3JV, etc. That is, all fractional occupations 
(plateaus corresponding to heights that are not integer 
multiples of the number of cavities N) disappear as k — > 
0, this compares nicely with the mean-field solution. The 
point at which this pinching occurs is called the critical 
chemical potential \i c . We find 
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[i c (n) = u + x{n) - x(n + I), 



(18) 
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FIG. 4: (Color online) These plots show the expectation 
value of the total number of excitations (L) as a function 
of (^ — uj)/P and n/j3, for two, three, four and five cavities 
in periodic boundary conditions, with A = 0. Note that the 
top boundary in each plot is the limit of calculations. The 
Hamiltonian matrix used to create each plot is truncated at 
/ = 12, 9, 8 and 5, respectively. 



ties, and independent of whether or not the mean- field 
approximation is used, as expected from the k — > limit. 

Although in general, one cannot analytically deter- 
mine the positions of all the boundaries, we find that 
for all topologies, the first boundary (between (L) = 
and (L) = 1) is described by the analytic equation 



fl — u> 



A + ZK 




A-; 



(19) 



where x( n ) 1S defined in Eq. and x(0) = — A/2. This 
is independent of the number or arrangement of cavi- 



where z is the number of nearest neighbors (the third 
column in Table [TJ) . Equation (|19p was determined by 
equating the smallest eigenvalue of with zero (the 
only eigenvalue of 7^°)). This has been compared nu- 
merically for z = 1, 2, 3, 4, 6 and is in excellent agreement. 
One cannot expect generic z analytic boundaries between 
(L) = 1 and (L) = 2 or higher to exist, indeed none 
were determined. This is because they are truly in the 
realm of many-body physics (unlike the lower boundary) . 
A higher boundary would be the solution of a commcn- 
surately higher order polynomial. In the simplest case, 
N = 2 and 1 = 2; from Eq. ([5]), s = 8, so an eighth order 
polynomial must be solved to determine the boundaries. 

For each plateau, the ground eigenstate can be calcu- 
lated (hence below we use "ground eigenstates" to refer 
to the different ground states for the different plateaus). 
These ground eigenstates are represented in the dressed 
state basis, even though many of the calculations above 
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FIG. 5: (Color online) The expectation value of the total num- 
ber of excitations (L) of two cavities with periodic boundary 
conditions in the ground state. The eigenstates of the four 
lowest plateaus are marked, (see Table JJ. The upper bound- 
ary marks the limit of calculations. 

were done in the bare basis. 

To easily represent this information, it is useful to de- 
fine two operators: the translation operator T (which 
shifts states to the right and moves the last state back to 
the beginning) and the permutation operator P m (which 
is the sum over T l applied to the state m times) 



states for some other band m, one needs to simply re- 
place |— , 1) by | — ,7n+l), and \g,0) by \ — ,m), in every 
instance. 

The differences between topologies are apparent when 
comparing, for example, the band m = 0, for 2 < (L) < 
N — 1, two excitations in a square topology compared 
with two excitations in a tetrahedron topology. In the 
square case, there is a different coefficient for excitations 
adjacent as to excitations separated, compared with the 
tetrahedron case, where all terms have the same coeffi- 
cient. 

Three different topologies (pentagon, six cavities with 
z = 3, and six cavities with z = 4) have some coefficients 
displayed to three decimal places, these have been cal- 
culated to twelve decimal places. An exact form is not 
derived, as these numbers represent solutions of polyno- 
mials of order > 50 no exact form is necessarily expected. 

We also examine the expectation value of the number 
of excitations of each cavity, (Li) and (L?), and find that, 
independent of topology, 

(U) = { 4 

so that the variance of Li is 



T(\si,m) 



|sjv-i,«iv-i) ® \sN,n N )) 



\s N ,n N ) ® |si,?7i) 



\sN-i,n N -i), 



P m (|si,ni) 



\sn, n N )) 



7H 



\sN,n N )) 



(20) 



(21) 



The ground eigenstates for all eleven topologies, up to 
(L) = N, i.e., total number of excitations equal to the 
number of cavities, are displayed in TableU See Fig.[5]for 
a diagrammatical example of how, for two cavities, the 
information from the table matches a phase diagram. 

Consider bands labeled by m, where 



iN < (L) < (to+ 1)N - 1 Vm = 0,l,. 



(22) 



One finds that the physics of each band has some striking 
similarities, hence the introduction of this parameter. 

In Table [I] only results for m = are displayed. 
While higher bands include many more possible states 
(e.g., |— , 1) ® |+, 1) <g> |™, 2)), we find higher bands have 
a surprisingly simple structure. To obtain the general 



var(ii) = v/(if)-(ii> 2 



(L) - mN ( (L) - mN 



(24) 



N 



N 



where the band 777 is defined by Eq. (|22|) . and both 
Eqs. ([23]) and are valid for i = 1,...,N. From 
this, one can determine that if (L) is an integer mul- 
tiple of N, var(Li) = 0, as expected. Also, if we con- 
sider the thermodynamic limit, where both the number 
of cavities and the number of excitations approach infin- 
ity (N — > 00, (L) — > 00), while the excitation density re- 
mains constant at p = (L)/N, we find that var(L^) — > 0. 

While this paper primarily focuses on phase changes 
as a function of k, one can also examine phase changes 
as a function of detuning A [T3|. Indeed, experimentally 
shifts in A may prove to be more accessible (via the Stark 
shift), as k cannot be changed post-fabrication in many 
systems. In Fig. [(^a) we plot (L) as a function of A//3 
and (/x — w)//3 for k = 0, and in Fig.[6]Jb) we do the same 
for k = lO^ 1 / 2 /?. There are fewer plateaus in (a) than 
there are in (b) because there are fewer plateaus (due to 
the pinching effect, discussed above) at k/(3 = 0. Note 
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FIG. 6: (Color online) These plots show (L) as a function 
of (/x — uj) /P and A//3 for two cavities in periodic boundary 
conditions with (a) k = and (b) k = 10~ 1//2 /3. The upper 
boundary in both cases marks the limit of calculations. Note 
that only even plateaus are present in (a); this is because of 
the pinching effect as n — » - in this limit, plateaus corre- 
sponding to fractional occupation do not exist. The white 
dashed line marks A//3 = 0, and aids the eye in seeing that 
the boundaries above the first (second) are symmetric (asym- 
metric) in plot (a) [(b)]. 



the symmetry around A = 0, in the second and subse- 
quent boundaries of Fig.[6ja) (c.f. Fig. 3 of 0), and the 
corresponding asymmetry in the third and subsequent 
boundaries of Fig. [B^b). This symmetry is perfect at 
k — 0, and the asymmetry increases with increasing k. 
A mean-field phase diagram that is comparable with 



the phase diagrams of the previous section can be made 
Q . An accurate comparison between the exact results of 
the previous section, with mean-field, is made when we 
consider topologies with z nearest neighbors with mean- 
field results for z nearest neighbors. In all eleven dis- 
tinct topologies tested, a very accurate match is seen. 
One such result is displayed in Fig. [7j We find that the 
boundaries from Eq. (|18p arc preserved in the mean-field 
solutions. 




i°g 1o 0</p) 



FIG. 7: (Color online) This plot shows the mean-field result 
for z = 1, overlaid with the exact cavity results for two cavities 
with one connection (white dashed lines), as in the fifth row of 
Table [T] Note the excellent agreement between exact results 
and mean-field approximation. 



In mean-field, the region with -0 = corresponds to the 
various Mott insulating lobes (e.g., \g,0), \— , 1), |— ,2) 
etc.), while ip > is the superfiuid state. The bottom 
lobe is described as the zeroth lobe (\g, 0)), the next lobe 
up as the first lobe (|— , 1)), and so on. In Fig. |5[a), we 
examine the underside of the first lobe in mean-field with 
z = 2. and overlay the boundary between (L) = N — 1 
and (L) = N for N = 2, 3, 4, 5, and 6 cavities in periodic 
boundary conditions. In Fig. [5Jb), we examine the un- 
derside of the second lobe in mean- field with z = 2, and 
plot the boundary between (L) = 2N — 1 and (L) = 2N 
for N = 2, 3, 4, 5 cavities, also in periodic boundary con- 
ditions. One can clearly see how as the number of cavi- 
ties increases, the boundaries approach the boundary of 
mean-field, and eventually may pinch off for each lobe 
entirely (so that plateaus of height (L)/N = 1, 2, . . . do 
not continue as k — > oo, but rather have finite size in this 
direction). These boundaries accord well with the struc- 
tures observed by Rosario and Fazio (Fig. 3 or Ref. [lol]). 
which were obtained independently by the density matrix 
renormalization group procedure, which lends weight to 
both quantum treatments. Furthermore, as N — > oo, our 
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exact results approach that of the mean-field, which have 
a more rounded cutoff for the Mott lobes than these finite 
cavity results. 




-3 -2.5 -2 -1.5 

iog 10 (K/p) 



-0.5 



-0.36 




-2.6 -2.4 
log 10 (K/p) 



FIG. 8: The hatched region represents the ip = area of the 
mean-field approximation for two nearest neighbors. The lines 
(or markers, for distinction purposes only) represent the exact 
calculations, (a) shows the boundary between the zeroth and 
first lobe, and (b) shows the boundary between the first and 
second lobe. Note that as the number of cavities increases, 
the lines tend to hug the upper mean-field lobe. This indicates 
qualitatively that as N — » oo, the exact calculations should 
approach the mean-field (TV = oo) limit. 



IV. DISORDER AND EFFECTIVE MODEL 
TEMPERATURE 

In this section, we first consider the modification of 
the chemical potential with small temperature increase 
(less than the scale for photon generation, kT <C Huj) 
and hence this modifies the phase diagrams above. We 
then examine fabrication disorder in the form of a nor- 
mal distribution of photon energies for each cavity. We 
show that this fabrication disorder is qualitatively similar 



to effective temperature, providing a connection between 
disorder and an effective temperature in this analogue 
system. 

Note that we set Boltzmann's constant ks = 1. We 
begin by differentiating the free energy F = E — TS with 
respect to the total number of excitations I, 



dF_ _ dE_ _ T dS_ _ dT 



(25) 



recalling the definition of chemical potential in Eq. (J7J), 
and assuming that temperature does not depend on the 
number of excitations (i.e., dT/dl = 0, assuming that 
the temperature scale is too low to generate a photon, 
i.e., kT <§C Huj), we get 



dE 



OS 
T — 

di ' 



(26) 



this then gives an effective chemical potential 



// = fJL + T 



dS_ 
dl ' 



(27) 



We calculate the entropy S in the n — > limit in the 
following manner. Assume that the photon blockade is 
complete, i.e., |— , 2, g, 0) -** | — , 1, — , 1), then we can con- 
sider each band [recall Eq. f22])] separately. More specif- 
ically, each band acts like a paramagnet [41j . Recall that 
a one-dimensional paramagnet is a line of spin states, 
where each spin can point up or down. Compare with 
our system, where each cavity can have either \g, 0), or 
|— , 1) (for band m = 0). Strictly speaking, each state 
(as in Table [T| is a pure state, and as such the entropy 
is zero. However, if we assume that the number of cav- 
ities N is very large, then the superposition states acts 
like a mixed state, and we can define entropy as for a 
paramagnet (essentially the logarithm of the number of 
microstates) by 



S(l) = In 



N 
I - mN 



(28) 



Note that this solution is only valid within each band, as 
such we can ignore the infinities that arise in dS/dl when 
I is an integer multiple of TV, as at these points the para- 
magnetic approximation does not apply. Recall that the 
phase diagram of (L) is concerned with finding the slope 
of the smallest energy eigenvalue with respect to fi, and 
that this slope is always an integer. The Hamiltonian of 
Eq. ([5]) is block diagonal; we know from earlier analy- 
sis that the ground state energy of each block Hy> has 
constant slope with respect to fj, of —I. Consider Fig. [31 
when temperature is included, each line will move, with 
respect to (// — oj)//3 by some amount to the left or to the 
right. For each /x, we choose the smallest energy eigen- 
value at that point, and take the negative slope at that 
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FIG. 9: (a) shows how the boundaries between plateaus 
change (when ft = 0) for 10, 100, and 1000 cavities with 
increasing temperature in natural units, (b) shows how the 
boundaries between plateaus change (when k, = 0) for two, 
three, and four cavities with increasing disorder, measured in 
units of standard deviation Note that when c = 0.4, the 
difference between the top line of the bottom pinch, and the 
bottom line of the top pinch, is (p — u))/ f3 =0.434, 0.361, and 
0.320 for two, three and four cavities, respectively. This gives 
some indication of a threshold of tolerance. 



point. For small finite temperatures, this is manifest as 
a "splitting" of the pinches, as seen in 0. We plot this 
splitting between the zeroth and first lobes, and the first 
and second lobe, for N = 10, 100, 1000 in Fig. g[a). 

Fabrication of a system of photonic cavities will un- 
doubtedly be subject to certain errors. Here we model 
uncertainty in the cavity frequency to. We assume that 
each cavity may be tuned individually to A,; = Vi = 
1,...,N (probably via the Stark shift), and as such 
model the Hamiltonian by 



N 

1=1 (29) 
+ Pi (erf a, + o-r a f) - re ^ a\ aj , 

(id) 

where the set {5\, S2, ■ • • , <5jv} is chosen from a normal 
distribution with zero mean and standard deviation <;. 
For fixed number of cavities and fixed <;, wc calculate 
1000 sets each of boundaries above I — to below I = 2N, 
and take the mean of the results. Results are shown in 

Fig. rib). 

One can see by comparing Figs, da) and db) that 
disorder and temperature produce qualitatively the same 
results. However, these disorder effects can only be calcu- 
lated up to four cavities due to the limitations of comput- 
ing resources, and the temperature analysis is only valid 
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FIG. 10: How T* , the temperature at which the top boundary 
of the first group in Fig.|9ja) meets the boundary of the second 
group, changes as a function of the number of cavities. 



for large numbers of cavities. Hence the two techniques 
cannot be compared directly. If the exact diagonaliza- 
tion technique could be extended to a larger number of 
cavities, it could be compared quantitatively with disor- 
der, matching properly the standard deviation with the 
effective temperature T. 

If one envisions Fig.^a) as temperature increases even 
further, there will be some temperature T* such that 
the top line from the bottom group (corresponding to 
I = N — 2) will meet the bottom line of the top group 
(corresponding to I = N + 1). Wc examine T* as a 
function of the number of cavities, and find that this 
is given by 



T* = (2 - V2) T(7V - 2) + T(N + 2) 
-T(N- 1) - F(N + 1) 



(30) 



where T(l') = I' ^|-L_,,, this function is plotted in 
Fig. 1101 T* appears to converge to a constant, non-zero 
temperature as N — > 00. 



V. CONCLUSIONS 

In this paper we present an intensive analysis of the 
Jaynes-Cummings-Hubbard model using the exact di- 
agonalization technique - studying the phase diagrams 
via the expectation value of the total number of excita- 
tions. We examine various topologies of small networks 
of cavities, and compare this work with the infinite cav- 
ity mean-field approximation, we find good agreement 



11 



in all topologies. We study the effective model temper- 
ature, and compare this qualitatively with disorder in 
the photon energy of the exact JCH, and also find good 
agreement. 
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